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A novel mathematical treatment is proposed for computing the time evolution of dynamic nuclear 
polarization processes in the low temperature thermal mixing regime. Without assuming any a 
priori analytical form for the electron polarization, our approach provides a quantitative picture 
of the steady state that agrees with the well known Borghini prediction based on thermodynamic 
arguments, as long as the electrons- nuclei transition rates are fast compared to the other relevant 
time scales. Substantially different final polarization levels are achieved instead when the latter 
assumption is relaxed in the presence of a nuclear leakage term, even though very weak, suggesting 
a possible explanation for the deviation between the measured steady state polarizations and the 
Borghini prediction. The proposed methodology also allows to calculate nuclear polarization and 
relaxation times, once the electrons/nuclei concentration ratio and the typical rates of the micro- 
scopic processes involving the two spin species are specified. Numerical results are shown to account 
for the manifold dynamical behaviours of typical DNP samples. 



I. INTRODUCTION 

Dynamic Nuclear Polarization (DNP) techniques are 
attracting increasing interest due to their proven ability 
of improving the sensitivity of Nuclear Magnetic Res- 
onance (NMR) experiments by several orders of magni- 
tude, not only in the solid state pQ but also in solution [2]. 
Although the disruptive potential of DNP is nowadays 
well known and accepted, the articulated physical sce- 
nario underlying these complex phenomena is still puz- 
zling the scientific community. 

In the solid state, three different DNP regimes can be 
specified according to the resonance frequency uo n of the 
nuclei to be polarized, the width Auj e of the Electron 
Spin Resonance (ESR) line of the free radicals used as 
polarization source and the magnitude of dipolar cou- 
pling between those radicals 1/T2 e . 

When uo n > Auj e the main mechanism for nuclei to 
polarize, known as Solid Effect (SE), proceeds via mi- 
crowave (MW) assisted forbidden transitions involving 
simultaneous flip-flops or flip-flips of one electron and one 
nucleus in mutual dipolar interaction. The SE model has 
been extensively described in [3 and the role played by 
spin diffusion was discussed in [4]. More recently, an ex- 
act quantum mechanical treatment of the SE based on 
the density matrix approach has been proposed by Hovav 
et al. [5]E] and by Karabanov and collaborators [7]. 

On the other side, when uj n < Au; e , nuclei can flip 
between different Zeeman levels within an energy con- 
serving three particle mechanism involving a simulta- 
neous flip-flop of two electron spins (referred as ISS 
process hereafter). Such process, driven by electron- 
electron and electron-nucleus time dependent dipolar in- 
teractions, does not involve forbidden transitions and it 
is thus in general more effective than SE. As far as the 
typical interaction time T2 e between different paramag- 
netic centres is long compared to the electron spin lattice 
relaxation time T\ e {i.e. the unpaired electrons are - on 



average - relatively far from each other), the polarization 
mechanism is referred as Cross Effect (CE) and was first 
analyzed in [5HTT j. The CE model has been successfully 
exploited to describe DNP with bi-radicals p~2] [13 , i.e. 
polarization procedures where the polarizing agents are 
tailored molecules carrying two unpaired electrons having 
resonance frequencies differing exactly by uo n . Similarly 
to SE, also the CE model has been rigorously computed 
by ab initio quantum mechanical techniques [14] [15] . 

However, the novel applications which in the last 
decade have renewed the attention on DNP due to their 
potential impact on biosciences [16] [17], are based on 
samples which generally fall in none of the two mod- 
els previously described. The nuclei of interest for this 
kind of applications (typically 13 C or 15 N) have a low- 
gyromagnetic ratio j n and, as a consequence, a resonance 
frequency uo n significantly smaller than the width of the 
typical radicals used as polarizing agent (trityls or ni- 
troxides). Moreover, the concentration of these radicals, 
usually above 10 mM in the solution to polarize, leads 
to non negligible electron spin-spin interactions, i.e. to 
the condition \jT<i e » that define the Thermal 

Mixing (TM) regime. Under TM assumptions, when a 
transition is saturated by an external, frequency selec- 
tive, microwave field, the whole electron spin distribu- 
tion reacts through energy conserving flip-flops (spectral 
diffusion) and evolves towards a new steady state dif- 
ferent from Boltzman equilibrium. Hence, via the pre- 
viously introduced ISS processes, the nuclear Zeeman 
populations are also perturbed and, depending on the 
frequency of the saturated electron transition, possibly 
result in depleting or enhancing the population of the 
ground state and thus in an enhanced nuclear magnetic 
order. At high temperature, where a linear expansion 
of the density matrix can be used, the evolution of the 
electron and nuclear polarizations is accounted by a set 
of rate equations formulated by Provotorov [3] [T8H2T] . 
At the typical temperatures (about 1 K) where most of 
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the DNP experiments aimed to obtain substantial nu- 
clear polarization values (> 0.1) are performed, however, 
a linear expansion of the density matrix is not allowed 
and consequently the Provotorov approach does not ap- 
ply. Fortunately in these conditions, the electron reso- 
nance lines of the free radicals used as polarizing agents 
are normally inhomogeneously broadened and can likely 
be depicted as the convolution of several individual pack- 
ets of given resonance frequency, mutually connected by 
the electron-electron dipolar interaction [3] . Hyperpolar- 
ization via TM at low temperature has been first dis- 
cussed by Borghini, who also calculated the steady state 
solution for the electron and nuclear polarization after 
imposing certain constrains to the model [22]. The same 
- steady state - result has been achieved by Abragam 
and Goldman [3 using a slightly different mathematical 
procedure. The Borghini prediction has the merit of re- 
producing, at least qualitatively, many observations ob- 
tained in DNP experiments at low temperature [23H25] . 
It has, at the same time, several limitations, the major 
of them listed below. 

i) It provides a picture of the steady state polariza- 
tion without contributing in any way to the under- 
standing of the time evolution of nuclear spin order 
towards equilibrium. 

ii) It is obtained by assuming an a priori analytical 
form for the electron spin polarization. 

Hi) It is derived under strong saturation conditions and 
in the limit of perfect contact between the electron 
and the nuclear reservoirs; the latter hypothesis 
consists in assuming highly effective ISS processes. 

As a consequence of the constrains ii) and iii) , the model 
leads to a final steady state nuclear polarization which 
substantially depends only on the lattice temperature, 
the magnetic field in which the DNP phenomenon takes 
place and the width of the electron resonance line. The 
sample specific parameters like electron and nuclear spin- 
lattice relaxation rates and the relative concentration of 
the two spin species play only a secondary role. 

In this work we aim to overcome these three limitations 
by introducing a dynamic analysis of the low temperature 

) based on a 
which spon- 



TM model (described in detail in Section [II 
set of rate equations presented in Section |III 
taneously provide the time evolution laws of the electron 
and nuclear polarizations without any a priori assump- 
tion on their functional form. The results obtained by 
numerically solving the rate equations are reported in 
Section IV and discussed in Section [V| with particular 
emphasis on the modifications occurring in the steady 
state and in the dynamic parameters in the presence of 
nuclear leakage and finite electron- nucleus exchange. The 
technical arguments underlying the derivation of the rate 
equations set are given in Appendix [A] and [Bj while in 
Appendix[Cj for convenience of the reader, the derivation 
of the Borghini equation for the steady state as proposed 
by Abragam and Goldman [3 is briefly recalled. 



II. MODEL DESCRIPTION 

A system made up of N n nuclear spins I with Larmor 
frequency uo n and N e electron spins S with mean Larmor 
frequency uo e is considered (S = / = 1/2). Typically uo e 
is about three orders of magnitude higher than uj n . Since 
all terms in the nuclear Hamiltonian with the exception 
of the Zeeman one are small in comparison to the other 
energy scales considered, nuclei are assumed to resonate 
all at the same frequency. 

The first assumption of the model is that the main 
contribution to the ESR line shape is the spread of g- 
factors (inhomogeneous broadening of the I type), that 
reflects the different orientation of the unpaired magnetic 
moments with respect to the external magnetic field due 
to single ion anisotropy. The practical relevance of this 
assumption is confirmed, for instance, by the experimen- 
tal evidence obtained from the widely used trityl radical 
family [23] . Inhomogeneously broadened lines are conve- 
niently decomposed in a sequence of narrow individual 
spin packets of frequency uji = uj e — A^, width 5uo and 
relative weig ht fi (see FIG. [TJ such that: 
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(i) 



valid for any ESR line. The continuous limit is recovered 
when Suj —> 0. For each electron packet, a local polariza- 
tion P ei = 2 (S l z ) may be defined, where () stands both 
for the quantum mechanic expectation value and the av- 
erage over all electron spins belonging to the i-th packet. 
The second assumption of the model is that uj n is smaller 
than the ESR line width, as actually is the case for low j n 
nuclei with almost any type of radical or for high gamma 
nuclei, such as 1 H, with broad line radicals (e.g. nitrox- 
ides). In view of the mathematical description proposed 
in the next section, it is useful to introduce at this point 
the variable 



Sn ri 



Suj ' 



(2) 



i.e. the number of electron packets corresponding to the 
nuclear Larmor frequency (see FIG.JT]). 

The third assumption is that the spin dynamics of the 
system is governed only by the five processes depicted in 
FIG. [2] and briefly explained here below. 

Microwave irradiation. Single electron transitions are 
stimulated by an external MW field at ujq = uj e — Aq that 
leads the electron spin system out of equilibrium (panel 
A of FIG. [2]). The characteristic time of this process is 
named Timw- For high microwave power Timw —> and 
the packet Aq is saturated. 

Spectral diffusion. The basic transition of this process, 
which involves only electrons and conserves both the total 
energy and the total electron polarization, is represented 
in FIG.[2j panel B. Spectral diffusion transitions are pro- 
moted by the dipolar interaction among electrons and 
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FIG. 1: Discretization of the ESR line. Each basic packet 
centered at uji — uo e — A* is characterized by a width 5uj and 
a weight f\. By way of example the irradiation frequency 
is set to ujq — uj e — Ao and the nuclear Larmor frequency 
corresponds to three basic packets (5n p = 3). 



are characterized by a time constant T2 e which in the 
thermal mixing regime is assumed to be much shorter 
than any other relevant time scale (in practice, in a solid 
solution used for DNP, it is typically < 1 fis [23]). 

ISS process. The mechanism is sketched in FIG. |2j 
panel C, and consists in a simultaneous flip-flop of two 
electron spins, belonging to packets separated by uo ni 
compensated by a nuclear spin flip. This process allows 
an energy transfer between nuclei and electrons and while 
the total energy and the electron polarization are con- 
served, the nuclear polarization is not. In a real system 
the ISS process, being promoted by dipolar interactions, 
involves only nuclei which are sufficiently close to para- 
magnetic centres. Thanks to nuclear spin diffusion, the 
local information is then spread throughout the entire 
system. A model which takes into account both these 
aspects (local dipolar interaction and spatial diffusion) 
would be very cumbersome to treat, because of the N n 
degrees of freedom. In the present work we adopt a mean 
field approach where all possible terns comprising a nu- 
cleus and two electrons flip with a characteristic time Tfgg 
which does not depend on the mutual distances between 
the three particles. Setting low values of Tiss the model 
mimics the fast spin diffusion limit, while setting high 
values of Tiss the model mimics the slow spin diffusion 
limit. 

Electron spin-lattice relaxation. The process, shown 
in FIG. [2] panel D, accounts for the contact between the 
electron system and the lattice which leads P e { (t) towards 



the electron thermal equilibrium polarization Pq: 

„ hu e . . 

P ° = - tanh 2^f' (3) 

where T is the temperature of the lattice. Its character- 
istic time constant T\ e ranges from hundreds of millisec- 
onds to few seconds [231 [26] , for T < 10 K. 

Nuclear spin-lattice relaxation ("leakage"). The pro- 
cess is shown in FIG. [2j panel E. The nuclear system 
is directly in contact with the lattice via slow processes 
(where electrons are not involved), with a characteristic 
time Ti n > 10 3 s, also called "leakage" terms. Analo- 
gously to the latter case, these processes lead nuclear po- 
larization, P n (t), towards the thermal equilibrium value 
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FIG. 2: Microscopic events in the thermal mixing model. 
Electron spins belonging to a generic packet i are represented 
in red with a subscript e, i, nuclear spins in blue with a sub- 
script n. Panel A: electron flip due to MW irradiation; Panel 
B: spectral diffusion (5 is a generic number of packets); Panel 
C: ISS processes; Panel D: electron spin- lattice relaxation; 
Panel E: nuclear spin-lattice relaxation ("leakage"), under the 
assumption Po,n ~ (see Eq.Q). 
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In the next section, the temporal evolution of the nu- 
clear polarization for the model described above will be 
determined by a closed set of equations involving only 
P n (t) and P e 4t). 



The numerical procedure can be summarized as fol- 
lows: 

PIS) gTT (Pe,i(t) LTE (P2A* + *t) 



Pn(t) 



Pn(t) 



A. STT 



Pn(t + dt) 



III. RATE EQUATIONS APPROACH 

To investigate the dynamic evolution of the model sum- 
marized in the previous section, we introduce a system 
of rate equations and solve them numerically. The dis- 
crete time step dt for the numerical integration is defined 
as the inverse of the sum of the single rates for the five 
processes sketched in FIG. [2] 

Microwave rate. The rate of all possible microwave 
events is W M w ^e/o/^iMW. 

ISS process rate. In the mean field approximation, the 
number of the possible processes involving two electrons 
and one nucleus is given by J2i fi N efi+5n p N e N n . The 
total ISS rate needs to linearly scale with the size of the 
system in order to assure a correct thermodynamic limit. 
To achieve this, as is usually done for fully connected 
models, the effective time constant of each three particle 
process must depend on the system size and scale as: 



r eff 
[ ISS 



Tissue N n , 



(5) 



where the constant Tjss is size independent. Then 
the total rate of all ISS events writes Wiss = 

Ysi fifi+5n p /T\ss. 

Electron and nuclear spin-lattice relaxation rate. The 
rate of all the electron spin-lattice relaxation processes 
is W e = N e /T\ e . Similarly the rate of all the nuclear 
spin-lattice relaxation processes is W n = N n /Ti n . 

Spectral diffusion rate. Being proportional to 1/T2 e , it 
is assumed to be by far the highest total transition rate 
among those characterizing the five considered mecha- 
nisms. A large rate makes the time step dt very small and 
consequently produces a dramatic slow-down of the sim- 
ulation procedure. To avoid that, the dynamic problem 
has been separated in two steps: a "Short-Term Ther- 
malization" (STT) step involving spectral diffusion and a 
"Long-Term Evolution" (LTE) where the system evolves 
according to defined differential equations under the ac- 
tion of the remaining processes. Depending on the value 
of Timw? the microwave pumping can be considered ei- 
ther as a fast process (contributing to STT) or as a player 
in the LTE. In the first case the time step is defined as: 



dt = 



W e + W lss + W n ' 
while in the second case one has: 

dt = 1 



W e + Wiss + W n + W 1M W 



(6) 



(7) 



The goal of this step is to "thermalize" the profile of 
the electron polarization P^i) at a generic time t, under 
the action of those processes considered as "fast". The 
steady state profile P e ,i(t) can be obtained imposing the 
detailed balance condition (see Appendix [A| : 



P Cj< (t) = -tanh [/3(Ai-c)]. 



(8) 



This equation depends on two parameters that can be 
computed using conservation principles. In this respect 
we need to discuss separately the two cases where spectral 
diffusion only or both spectral diffusion and microwave 
promoted processes are viewed as fast. 



Spectral diffusion 



This case is relevant for relaxation experiments (with 
microwave field off by definition) or when a non satu- 
rating microwave field is applied in a DNP experiment. 
To calculate the two parameters c and /?, the conserva- 
tion of both the total energy and the total polarization is 
imposed: 



J2fAi {Pe,i(t)-PlS))=0 

Y J fi{Pe,m-pls)) = Q- 

2. Spectral diffusion and saturating microwaves. 



(9) 
(10) 



For high irradiation power, the microwaves act as an 
infinite bath for energy exchange and force the polariza- 
tion P e ,o(t) to and in turn, by Eq.([8|, c = A . In order 
to evaluate the second parameter p7the output of the 
STT step is conveniently written as: 



Pe,0(*) = Pefi(t) + tPefl 



Hi ^ 
v5p mw 



(11) 
(12) 



where 5P e% i and 5P MW are the variation of the polariza- 
tion of the i-th packet induced by spectral diffusion and 
by microwaves respectively. The conservation of the total 
energy and of the total polarization can now be written: 



^2f i A i 6P e , i = 

J2fiSPe,i=0. 



(13) 
(14) 
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Eg.|l4| can be recast as f 5P ei0 = ~ S;/o f^ P e,i^ so 
that: 

= ^/,(A,-A )^Pe,z=0. (15) 

Summing on both sides fi(^i ~ Ao)P® i (t) we S e ^' 

X>(A< - A )P e ,i(t) = ]>>( A ' " A o)^ i(*)- (16) 
The condition for f3 

£/i +tanh[/3(A i - A )]} (A,-A ) = (17) 

identifies the unique solution for P e ^(t). 

B. LTT 

In this second step the profile of the electron polariza- 
tion P^^t + dt) at time t + dt is deduced from the output 
of the STT. When microwaves are off (relaxation) or they 
have been already taken into account in the STT, the rate 
equations for P Gyi and P n , as discussed in Appendix [Bj 
can be set down in the form: 



Pe,i(t + dt) = P e ,i(t)+ 



dt 



Tie 



2Piss 



(18) 



P n (t + dt) = P n {t) + 
Pon — Pn{t) 



dt 



2T lss N r 



fifi 



i-\-8n v ^-n 



where II_ = II_(z,t),II + = II+(z,t),n n = U n (i,t) are 
given by the expressions: 

n_ = P e ,i-6n p (t)-P ei i(t)+ P n(t) [l-Pe,i-6n p (t)P ei i(t)] 
n+ = Pe,i+5n p (t)-P e 4t)+P n (t) [l-Pe,i + 5n p (t)P e:i (t)] 

n n = P e ,i+6n p (t)-P et i(t)+P n (t) [l-P e ,i+Sn p (t)P e ,i (t)] . 

The effect of a non saturating MW irradiation can be 
included by the term — P e ,o(^)/TiMW i n the rate equa- 
tion for P e5 o(t), describing the electron polarization of 
the irradiated packet. 

IV. NUMERICAL RESULTS 

In this section we exploit the formalism described 
above to calculate, in specific situations, those curves 
which are normally measured in DNP experiments, 
namely build-up and relaxation curves. The so-called 
build-up curve describes the polarization growth over 



time under selective microwave irradiation and it is char- 
acterized by two parameters, the final nuclear polariza- 
tion P n : 



P n = P n (t oo), 



(19) 



and the characteristic time T po \ that we define by the 
equation: 



Pn(t = T vol )=P n [l 



1 



(20) 



with P n (t = 0) = 0. Tp i corresponds to the usual time 
constant in case of exponential build-up but, through 
Eq.(|20|), one has a more general quantification of the 
growing speed regardless the specific form of the polar- 
ization function. 

Relaxation curves, on the other hand, represent the 
spontaneous (at microwaves off) equilibrium recovery of 
a system that is prepared out of equilibrium at t — 0. The 
characteristic time is named T re i ax and, if one starts with 
a fully polarized system (P n (t = 0) = 1), the constant 
Prelax is defined by: 



Pn(t — Treiax) — Po,n ( 1 



1 



(21) 



where Po, n is defined in Eq.Q. 

Setting of the parameters. In order to make the re- 
sults of our simulation relevant for a deeper understand- 
ing of experimental data, a number of basic parameters 
are borrowed from some of the best known DNP sam- 
ples for biomedical applications. In this framework, the 
typical number of 13 C nuclei (that in most cases are the 
object of DNP) is in the range 5 - 15 M, while the radical 
concentration is 10 - 20 mM [27]. Is thus reasonable, for 
our purposes, to set N n /N e = 1000. Concerning the dis- 
tribution of the electron spin resonances, it looks worth 
to use a Gaussian function with a full width at half maxi- 
mum Auj e = 63 MHz, as surrogate of the ESR line of the 
trityl doped samples studied in [23] . This is actually one 
of the few ESR line widths which have been measured in 
a standard DNP environment (a magnetic field of 3.35 T 
and a temperature of 1.2 K) and that turned out to be 
inhomogeneously broadened by ^-factor spreading. Fi- 
nally, despite the fact that our approach is not restricted 
to the case of saturating microwaves, we work under the 
approximation T\mw ^ oo to better mimic what is nor- 
mally done in actual experiments. 

ESR line discretization. The Gaussian function used 
as model of the ESR line is truncated at 3<r, where <r = 27 
MHz. This defines a frequency interval of about 160 MHz 
to be covered by electron spin packets. Three different 
discretizations of the line are employed, with different 
packets widths: Sou — 32, 10.6 and 6.4 MHz. The cor- 
responding number of packets N p is thus equal to 5, 15 
and 25. The value of uo n is set at 32 MHz (resonance fre- 
quency of 13 C nuclei at the typical DNP magnetic field 
of 3.35 T) which corresponds respectively to Sn p = 1,3 
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FIG. 3: Polarization build-up curves at 3.35 T and 1.2 K 
of a DNP system characterized by the following parameters: 
N n /N e = 1000, Ti e = 1 s, Tiss = 0.001 s, T ln -> oo, for 
different discretizations of the ESR line: N p = 5 (blue line), 
N p = 15 (red line), N p = 25 (gray line). The final nuclear 
polarization is P n = 0.825 reached with a characteristic time 
constant T po i = 1140 s. 



and 5 for the three selected values of N p . The MW fre- 
quency is set equal to Uq = uo e — a, which corresponds to 
saturating the packet 2 (for N p = 5), 5 (for N p = 15) and 
8 (for N p — 25). For this value of c^o the final nuclear 
polarization is maximal. 

The results we show are obtained by integrating the 
rate equations system (18) together with the STT step 
defined by the detailed balance condition (J8|, with c = 
Aq and j3 set by Eq.(17). The boundary conditions of 



FIG. 4: Relaxation curves at 3.35 T and 1.2 K of a DNP 
system with N n /N e 1000, T le 1 s, Tiss 0.001 s, T ln -> 
oo, when N p = 5 (blue line), N p = 15 (red line), N p = 25 
(gray line). The relaxation time is about 15600 s. 



curve with N p = 5 is already representative of the con- 
tinuous limit. The characteristic time T re i ax « 15600 s, 
evaluated by means of Eq.(21) with Po, n = 0, turns out 
to be about 10 times longer than the corresponding T po \. 



the ESR line shape are correctly implemented imposing 

fi+5n p = (fi-6n p = 0) if 2 + Sfl p > Np (l - Sfl p < 0). 



In FIG. [5] the effect of nuclear leakage is investigated. 
The thermal contact between the nuclear spin system 
and the lattice significantly affects the build-up curve 
even when T in is much slower than any other transition 
considered. For T\ n = 10000 s, the final nuclear polar- 
ization goes down to P n = 0.707, which corresponds to a 
reduction of 15% with respect to P n in absence of leak- 
age, while the corresponding polarization time becomes 
T pol = 910 s. 



A. Fast ISS limit 

The case of highly effective contact between electrons 
and nuclei, which corresponds to very fast ISS processes 
compared to electron spin lattice relaxation (Tiss <^ Ti e ), 
is considered first, by setting Tiss = 10 _3 Ti e and T ie = 1 
s (the latter from [23 ) in our computational tool. 

In FIG. [3j three different build -up curves, computed 
with an increasing value of N p and in absence of leak- 
age (Ti n —> oo) are presented. The continuum limit is 
approached very fast and a good convergence is reached 
already with N p = 5. The final nuclear polarization is 
P n = 0.825 and the characteristic time constant of the 
polarization build-up curves T po i is equal to 1140 s. 

In FIG.|4j three relaxation curves, describing the evo- 
lution of the nuclear polarization in the absence both of 
microwaves and leakage and calculated with a different 
discretization of the ESR line, are shown. Again, the 



B. Competing ISS and electron spin lattice regime 

The case of finite contact between electrons and nuclei, 
Tiss ~ Ti e , is now considered. 

In FIG. [6] and FIG. [7] build-up and relaxation curves, 
in absence of leakage (Ti n — >• oo) for T\ e = 1 s and Tiss 
within the range 0.001 - Is, are analyzed. Two regimes 
can be identified. For fast ISS process both T po i and 
T r eiax are not affected by the particular value of Tiss, 
since the bottleneck process is represented by the electron 
spin-lattice relaxation. Conversely, when Tiss grows the 
ISS process becomes the rate determining step and both 
polarization and relaxation times are enhanced. With 
Tiss = 0.1 s for instance, one obtains T po i = 5300 s and 
T r eiax = 44900 s. The ratio between the two character- 
istic times remains constant and close to 10. Finally no 
influence on the final polarization P n = 0.825 is observed 
on varying T^s- 

Surprisingly, when an almost negligible leakage term 
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is introduced (T in = 10000 s), the final nuclear po- 
larization becomes strongly dependent on the effective- 
ness of the contact between nuclear and electron systems 
(FIG. [8|. For Tiss = 0.1 s, the final nuclear polariza- 
tion is P n = 0.529, which corresponds to an important 
reduction of 36% that could be relevant for explaining 
some experimental observations. The polarization time 
is measured equal to 3200 s. 

A comprehensive summary of the role played by the 
different parameters is represented in FIG. [9] and FIG.flOl 
Two distinct regimes are clearly observable: highly and 
poorly effective contact between nuclear and electron sys- 
tems. The former (Tiss —> 0) is characterized by high 
levels of polarization and short build-up times, both de- 
pendent on Ti e , whilst the latter (Tiss oo) has low 
levels of polarization and long build-up times, both inde- 
pendent on Ti e . Clearly, in the first regime, the bottle- 
neck role is played by Ti e , while Tiss is sufficiently high 
not to affect the system dynamics. In the second regime 
their role is reversed. 



V. DISCUSSION AND CONCLUSION 

To date, the main attempt to give a theoretical descrip- 
tion of the DNP phenomenon at low temperature in the 
TM regime is the prediction proposed by Borghini [3l [22] 
for the steady state nuclear polarization, that makes use 
of the standard parameterization: 

P n = tanh (p n u n ) with/3 n = H/(2k B T n ), (22) 

where T n represents the temperature reached by the nu- 
clear reservoir under microwave irradiation. To compute 




6000 



FIG. 5: Polarization build-up curves at 3.35 T and 1.2 K of 
a DNP system with N n /N e = 1000, T le = 1 s, Tiss = 0.001 
s for different nuclear intrinsic relaxation times: T\ n — > oo 
(blue line), T ln = 10000 s (red line), T ln = 1000 s (gray line), 
Ti n = 100 s (orange line). 
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FIG. 6: Polarization build-up curves at 3.35 T and 1.2 K 
of a typical DNP system with N n /N e = 1000, T le = 1 s, 
Ti n — > oo for different values of contact time Tiss (0.001 s 
(solid blue line), 0.01 s (solid red line), 0.1 s (solid gray line), 
0.2 s (dashed blue line), 0.5 s (dashed red line), 1 s (dashed 
gray line)). 




200 000 



400 000 600 000 

t[s] 



800 000 



lxlO 6 



FIG. 7: Relaxation curves at 3.35 T and 1.2 K of a DNP 
system with N n /N e 1000, T le = 1 s, T ln -> oo for different 
values of contact time Tiss (0.001 s (solid blue line), 0.01 s 
(solid red line), 0.1 s (solid gray line), 0.2 s (dashed blue line), 
0.5 s (dashed red line), 1 s (dashed gray line)). 



this temperature, the original derivation conjectures that 
the final electron polarization profile (P e ,i = P e ,i(t 
oo)) takes the form: 

P e , i = -tanh[/3(A i -c)]. (23) 



Actually, when the spectral diffusion is the fastest pro- 
cess, this form can be rigorously derived - without any 
ad hoc conjecture - by imposing the detailed balance con- 
dition as shown in Appendix [Aj In addition to Eq.(23), 
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a perfect thermal contact between the nuclear and the 
electron system and an infinite microwave power were 
assumed in the approach proposed by Borghini, so that 
P = P n and c = Aq. 

Under these assumptions, the celebrated relation (re- 
derived in Appendix [C] for convenience of the reader): 



1.0 



£)/i(A<- Ao)P e ,i + A Po 



u n §^P n =0 (24) 



gives a unique solution for T n and thus for P n . In 
the absence of leakage, the value of P n predicted by 
Eq.(24) at given temperature and field depends only on 



the ESR line shape and the irradiation frequency. For a 
Gaussian shape the maximal enhancement corresponds 
to ujo — uj e — <j, where a is the standard deviation of 
the electron frequency distribution. The Borghini predic- 
tion is in good qualitative agreement with many exper- 
imental observations obtained in low temperature DNP 
experiments. However it does not provide an accurate 
quantitative description of experimental data and in par- 
ticular of the maximum final nuclear polarization that 
turns always out to be overestimated [23, 28-30 . When 
a typical trityl doped sample, with an EPR line width 
of 63MHz, Ti e = 1 s, and N n /N e = 1000 is considered, 
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FIG. 8: Polarization build-up curves at 3.35 T and 1.2 K 
of a DNP system with N n /N e = 1000 and T le = 1 s for 
different values of contact time Tiss (0.001 s (solid lines) and 
0.1 s (dashed lines)) and nuclear intrinsic relaxation times T\ n 
(— > oo (blue lines), 10000 s (red lines) and 1000 s (gray lines)). 



0.6 



0.4 



0.2 



0.0 



io- 



0.001 



0.01 
T,ss [s] 



0.1 



FIG. 9: Final polarization versus contact time Tiss for dif- 
ferent values of the electron relaxation time T\ e (2 s (circles), 
1 s (squares) and 100 ms (triangles)). Remaining parameters 
are set as follows: N n /N e = 1000, T ln = 10000 s. Dashed 
lines represent polarization final values as estimated by the 
Borghini prediction in Eq.(24). 



Eq.(24) predicts a maximum steady state polarization 
P n = 0.825 in absence of leakage. The same equilibrium 
value is obtained by our mathematical approach (FIG. [3] 
and FIG. [6]) under the same physical constrains. 

On the other hand the polarization levels achieved ex- 
perimentally are in general much lower (see for instance 
Ref. [23[ [27] and [31]). To justify this discrepancy, one 
can call on dissipation terms. In fact, if a finite T\ n is ^ 
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FIG. 10: Polarization build-up times as function of the con- 
tact time Tiss for different value of the electron relaxation 
time Ti e (2 s (circles), 1 s (squares) and 100 ms (triangles)). 
Remaining parameters are set as follows: N n /N e = 1000, 
Ti n 10000 s. 



considered in Eq.(24), a loss of P n is obtained, that how- 
ever, with Ti n generally being longer than 10 4 s [23], can 
hardly exceed 10 - 15% (cfr. FIG. [§. A severe reduc- 
tion of P n is obtained instead by means of the equation 
set (18), introducing a finite contact between nuclei and 



electrons in presence of a small leakage (see FIG. |8j). It 
is worth to notice that in absence of leakage, i.e. when 
nuclei are isolated by the lattice, the final polarization 
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is not affected by the efficiency of the contact between 
electron and nuclear systems. 

Besides providing a more flexible scenario for P n , the 
proposed mathematical framework allows the computa- 
tion of polarization and relaxation times. Two regimes 
have been identified. 

When the contact between nuclei and electrons is 
highly efficient, the bottleneck of the spin dynamics is T\ e 
and, by setting the parameters of the simulation accord- 
ing the experimental conditions used in [23] one obtains 
Ppoi ~ 10 3 s and T re i ax w 10 4 s, as effectively measured in 
[23] . It is interesting to observe that, in this fast exchange 
limit, an estimation of the order of magnitude of T po i 
and T re i ax can be heuristically derived as follows. During 
polarization, the nuclear system transfers energy to elec- 
trons, which are cooled by the lattice. In a time Ti e , the 
lattice can absorb an energy proportional at most to 7V e , 
so that T pol w N n T le /N e w 10 3 s, with N n /N e = 1000. 
During relaxation, the effective number of electrons that 
can absorb energy from the nuclear system, proportional 
to (1 - P e ,i)(l + P e ,i) ~ (1 - P§), collapses to « 0.1 due 
to the high electron polarization P e i w Pq = 0.95 at 3.35 
T and 1.2 K, thus explaining the factor T re i ax /T po i « 10 
observed both in simulations and experiments. 

On decreasing the effectiveness of the contact between 
nuclear and electron spins, a different regime is estab- 
lished, where the ISS process becomes the rate deter- 
mining step and both T po i and T re i ax become longer. This 
provides a possible explanation of why T po i can range be- 
tween very different values in samples with the same ratio 
N e /N n polarized in analogous conditions. By way of ex- 
ample, one can compare the [l- 13 C]-pyruvic acid samples 
doped with 15 mM of trityl radical studied in [23] [27] . 
where N e /N n = 1000 and T pol = 1200 s with the [1- 
13 C] -butyric acid sample mixed with 20% in volume of 
DMSO and doped with 10 mM trityl radical analyzed in 
[31], having a much longer T po i = 3400 s while the ratio 
N e /N n remains almost unchanged. 

In conclusion, we propose a novel approach based on 
rate equations for studying the dependency of dynamic 
nuclear polarization in the low temperature thermal mix- 
ing regime (T2 e <C T\ e ) from the microscopic transitions 
involving electron and nuclear spins. This approach al- 
lows the recovery of the whole build-up curve and, in 
the limit of perfect contact between nuclei and electrons 
and infinite microwave power, leads to the same final nu- 
clear polarization predicted by Borghini. In addition, by 
tuning the efficiency of the exchange interaction between 
nuclei and electrons, different values of P n are reached, 
providing an interpretation key for those experimental 
observations of P n which are not simply accounted for 
by leakage terms depending only by spin concentration 
and spin-lattice relaxation times. 

The rate equation approach can be easily extended to 
more complex experimental systems. A second nuclear 
reservoir which also participates in TM could e.g. be in- 
cluded to interpret the dynamic experimental data mea- 
sured in nitroxyl doped samples (where both 13 C and 



1 H Larmor frequencies do not exceed the ESR line width 
[24] [28] [29] ) or, similarly, in trityl doped samples contain- 
ing 13 C and 89 Y nuclei, both in contact with the electron 
reservoir [25] [30] . Finally the versatility of the approach 
proposed here would easily allow to introduce new dis- 
sipative processes violating the precise assumptions and 
conservation principles the Borghini prediction is based 
on and that could possibly be useful to justify the many 
unexplained observations of low temperature DNP, e.g. 
the reduction of P n on increasing N e [27] [32] . 



VI. ACKNOWLEDGEMENT 

Christophe Texier and Thomas Guedre are grate- 
fully acknowledged for the fruitful discussions. This 
study has been supported in part by Regione Piemonte 
(POR FESR 2007/2013, line 1.1.1), by the COST Ac- 
tion TD1103 (European Network for Hyperpolarization 
Physics and Methodology in NMR and MRI) and by 
ANR grant 09-BLAN-0097-02. 



Appendix A: Detailed balance 



After defining P~^ { as the fraction of electrons up and 
P~ i as the fraction of electrons down belonging to the 
packet z, the detailed balance condition under the process 
depicted in FIG. [2] panel B, writes: 



P(t,i-5 C^e,i) Pe,i+8 ~ ^e,i-S C^e[i) ^e,i+cV 



(Al) 



Then, by using the relation 



P+ = 



P~ = 



(1 + Pe,i) 

2 

(1 - Pe,i) 



one comes to an equation for the electron polarization 

(1 + P e>i _«) (1 - P e4 ) 2 (1 + P e>i+<5 ) = 



= (1 



5) (1 + P e ,i) (1 " Pe,i+S) , (A2) 



that can be solved in the continuum limit where Suj — >• 0, 
N p — >• oo and P e i — »• P e (A^). It is sufficient to con- 
sider the transitions between consecutive packets, so that 
P e ,i+s = Pe,i+i Pe(^i + and write a second order 
expansion 

P e (A, + Slj) « P e (A z ) + SuP^Ai) + — P e "(A,). (A3) 



By combining relations (A2) and fA3j), one gets the sec- 



ond order differential equation 

2P e (A i )P^(A i ) 2 + i^(Ai) (1 - P e (A i ) 2 ) = 0, (A4) 
whose general solution 

P e (Ai) = - tanh (/3(Ai - c)) (A5) 



10 



was used in our treatment as starting point for what we 
called "STT" process. The same parametric function has 
been postulated in the derivation of the Borghini predic- 
tion for the steady state. 



Appendix B: Microscopic derivation of rate 
equations 

The term proportional to 1/Tiss and describing the 
ISS process, used in the rate equations (18), is here de- 



rived for the electron polarization first and then for nu- 
clear polarization. 



1. Electron polarization 

Be the fraction of electrons up belonging to the 
packet z, P+ the fraction of nuclei up, P~ i the fraction 
of electrons down belonging to the packet z and P~ the 
fraction of nuclei down. When the ISS event depicted 
in FIG. [TT] occurs, the fraction of electrons up in the z-th 
packet is reduced by l/{N e fi). The number of possible 
transitions is the product of: 

• the number of the electrons up in the z-th packet: 

• the number of the electrons down in the (i+5n p )-th. 
packet: N e f i+S n p P~ i+5 n^ 

• the number of nuclei down: N n P~ . 

The rate of such process is l/(TissN e N n ), and the total 
reduction of Pj~- in the time interval dt is: 



dt 



(Bl) 



2iss' 

The total variation of P^ { induced by all possible ISS 
transitions, SP^, is given by: 

dt 



3P(t,i ~ jt fi+5n p (^P e ^P^i+5n p Pn ^i^i+Srip^n ) 
+ fi-5n p (P e ^P^i-S np P n ~ Pe~,iP e ,i-5n p Pn^ ' 

Using the relations: 



e,i 



(1 + Pe,i) 



(1--Pe,i) 



+ = (1 + Pn) = (1 ~ Pn) 

2 ' 2 ' 

the total variation of P e ^ induced by all possible ISS 
processes, oP e ,» = 2oP/-, can be written as follows: 



6P e ,i 



dt 



{f i+Snp [(l-P eii )(l+P eA+Snp )(l+P n ) 



4Tiss 

-(l+P e ,i)(l-P e , i+4np )(l-P n )] 
+fi- Snp [(1-P e>i )(l+P e>i _ 5 „ p )(l-P„) 

-(l+P e , i )(l-Pe,i-fc p )(l+Pn)] } • (B2) 
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e,i e,i+Sn p n 



e, i e, i+8n p n 
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FIG. 11: Schematic representation of one possible ISS event. 



The term proportional to 1/Tiss m the fir st e quation of 
set ( 18 ) can be now easily derived from Eq.(B2 ) by means 



of simple algebraic calculations. 



2. Nuclear polarization 

When the event depicted in FIG. [TT] occurs, the frac- 
tion of nuclei up is increased of a factor 1/N n . The num- 
ber of possible transitions is the product of: 

• the number of the electrons up in the z-th packet: 

NefiP+i, 

• the number of the electrons down in the (i + 5n p )-ih 
packet: N e f i+Snp P~i + Sn^ 

• the number of nuclei down: N n P~ . 

The rate of such process is l/(T\ssN e N n ), and the rele- 
vant increment of P+ in the time interval dt: 



N e dt 
N n T l3S 



fifi+5n p Pt,iP e ,i+5n p Pn ' 



(B3) 



Considering now all the possible processes, the total vari- 
ation of P+ induced by the ISS process, £P+, is given 
by: 



6P+ 



N e dt 
N n T lss 



^ ^ fifi+5n p 



p+ p~ p~ 

e,i e,i+5n p n 



_p- p+ p+ 

e,i e,i+5n p n 



(B4) 



Following the line previously described for P e ^, one im- 
mediately arrives to the equations for P n reported in the 
main text (18). 



Appendix C: Borghini relation 

To facilitate the reading of the manuscript, we report 
the derivation of the Borghini relation in Eq.([24]) accord- 
ing to the line proposed in [3]. The energy of the whole 
electron and nuclear system is conveniently split into two 
reservoirs: the Zeeman electron contribution 



EzM = l^eflLJeJ^fiPeAt) 



(CI) 
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and the non-Zeeman electron plus Zeeman nuclear term 

E N z-z n (t) = ^N n hw n P n (t) - ^N e hJ2 A ifi p eAt)- 

(C2) 

The time evolution of the two energy reservoirs is de- 
scribed by the following equations: 

dP e ,i(t) 



2 e ^ Jl dt 



dt 



(C3) 

' — ■ a 6 

dt 2 n n 2 ^ u * dt ' 

and since either the spectral diffusion and the ISS pro- 
cess conserve both Ez e (t) an d ENZ-z n (t), one obtains: 



dE Ze (t) l AT . 

= -N e HuJ e 

dt 2 



E,P - Pe,i{t) n Pe,0 
Ji ^ JO 



^NZ^ = 1 

dt 2 



^1MW 
^e,0 



. J-le J-1MW 



1 AT * P ^ 
^ J- In 



(C4) 



It is important to observe that the evolution of the two 
energy reservoirs depends on the time progression of all 
P e ,i(t) and P n (t), that is the full solution of the system 
of rate equations reported in Eq.([l8|. As far as only the 
steady state solution is required, however, it is sufficient 
to impose the simultaneous vanishing of both right-hand 
sides of Eq.(C4). Thus, by multiplying the first Eq.(C4) 
by Ao/c^e and adding it to the right-hand side of the 



second Eq.(C4), one gets rid of the microwave transition 
probability 1/Timw and obtains the celebrated Borghini 
relation given in Eq.(24). 
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